Onset of sliding friction in incommensurate systems 
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We study the dynamics of an incommensurate chain shding on a periodic lattice, modeled by 
the Frenkel Kontorova hamiltonian with initial kinetic energy, without damping and driving terms. 
We show that the onset of friction is due to a novel kind of dissipative parametric resonances, 
involving several resonant phonons which are driven by the (dissipationless) coupling of the center 
of mass motion to the phonons with wavevector related to the modulating potential. We establish 
quantitative estimates for their existence in finite systems and point out the analogy with the 
induction phenomenon in Fermi-Ulam-Pasta lattices. 
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The possibility of measuring friction at the atomic 
level provided by the Lateral Force Microscopes [Q and 
Quartz Crystal Microbalance |^ has stimulated intense 
research on this topic Phonon excitations are the 
dominant cause of friction in many cases Most stud- 
ies are carried out for one-dimensional non-linear lat- 
tices p[-p2| and in particular for the Frenkel-Kontorova 
(FK) model where the surface layer is modeled 

by a harmonic chain and the substrate is replaced by 
a rigid periodic modulation potential. The majority 
iP-p^ examines the steady state of the dynamical FK 
model in presence of dissipation representing the coupling 
of phonons to other, undescribed degrees of freedom. 

We study the dynamics of an undriven incommensu- 
rate FK chain. Our aim is to ascertain whether the ex- 
perimentally observed superlubricity can be due to 
the blocking of the phonon channels caused by an incom- 
mensurate contact of the two sliding surfaces. Therefore 
we do not include any explicit damping of the phonon 
modes, since we wish to find out if they can be excited at 
all by the motion of the center of mass (CM). In an earlier 
study, Shinjo and Hirano ||^ found a superlubric regime 
for this model, where the chain would slide indefinitely 
without dynamic friction but with a recurrent exchange 
of kinetic energy between CM and a single internal mode. 
We will show that their finding is oversimplified by either 
too short simulation times or too small system sizes. The 
inherent non-linear coupling of the CM to the phonons 
leads to an irreversible decay of the CM velocity, albeit 
with very long time scales in some windows. The dissi- 
pative mechanism is driven by the coupling of the CM 
to the modes with modulation wavevector q or its har- 
monics, LOnq^ and consists in a novel kind of parametric 
resonances with much wider windows of instabilities than 
those deriving from the standard Mathieu equation . 
The importance of harmonic resonances at Unq has been 
pointed out before , with the suggestion that 

they could be absent in finite systems due to the dis- 
creteness of the phonon spectrum. However, it has not 
been realized that they act as a driving term for the onset 



of dissipation via subsequent complex parametric excita- 
tions which we shall describe, establishing quantitative 
estimates for their existence in finite systems. A related 
mechanism has recently been identified in the resonant 
energy transfer in the induction phenomenon in Fermi- 
Ulam-Pasta lattices [ p^ . 

We start with the FK hamiltonian 
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where u„ are the lattice positions and I is the equilibrium 
spacing of the chain for A = 0, A being the strength of 
the couphng scaled to the elastic spring constant. We 
take an incommensurable ratio of I to the period m of 
the periodic potential, namely m = 1, ? = (1 -I- \/5)/2. 
We consider chains of N atoms with periodic boundary 
conditions un+i = Nl + ui. Hence, in the numerical 
implementation, we have to choose commensurate ap- 
proximations for I so that I x N = M x 1 with N and M 
integer, i.e. we express / as ratio of consecutive Fibonacci 
numbers. The ground state of this model displays the so- 
called Aubry transition from a modulated to a pinned 
configuration above a critical value Ac = 0.14. Here we 
just note that in the limit of weak coupling (A << Ac), 
deviations from equidistant spacing I in the ground state 
are modulated with the substrate modulation wavevec- 
tor q = 2TTl/m [ p^ as due to the frozen-in phonon ujq. 
Higher harmonics nq have amplitudes which scale with 
A". 

We define the CM position and velocity as Q = 
7fJ2n'^n,P=j;J2n Pn- By Writing u„ = nl + Xn + Q, 
the equations of motion for the deviations from a rigid 
displacement x„ read 

Xn — Xn+i + Xn-1 — + A COS {qn + 211x^1 + ZttQ) (2) 

We integrate by a Runge Kutta algorithm the N Eqs. (|^) 
with initial momenta p„ = Pq and Xn{t = 0) correspond- 
ing to the ground state. For a given velocity P, parti- 
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cles pass over maxima of the potential with frequency 
O = 2'kP, the so-called washboard frequency PJlQ]. 
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FIG. 1. Time dependence of the CM velocity for several 
values of Po for = 144, A = 0.05 ~ ^. The dashed line 
corresponds to Pq = ^ — 0.2966. Close to higher reso- 
nances (solid dots) a similar oscillatory behavior is observed, 
accompanied by a slower decay which is not apparent on the 
timescale of the figure. 



mentum for A ^ Ac/3 and several values of Pq- According 
to the phase diagram of Ref. a superlubric behavior 
should be observed for this value of A and Pq > 0.1. We 
find instead a non trivial time evolution with oscillations 
of varying period and amplitude and, remarkably, a very 
fast decay of the CM velocity for Pq ^ ^q/ (27r) despite 
the absence of a damping term in Eq. (^. A similar, but 
much slower, decay is found for uPq ~ w„q/(27r). In the 
study of the driven underdamped FK [|j it is shown that, 
at these superharmonic resonances, the differential mo- 
bility is extremely low. Here, we work out an analytical 
description in terms of the phonon spectrum which ex- 
plains this complex time evolution and identifies the dissi- 
pative mechanism which is triggered by these resonances. 
In the limit of weak coupling A it is convenient to go from 
real to reciprocal space by defining Fourier transformed 



coordinates Xk = jf J2n ^ 
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where k — 27rn/A^ and the normalization is chosen to 
remove the explicit A^-dependence in the equations of 
motion, which become: 

Xk = -ojlxk + ^Y. [e^«"e'2.Qg.2..„ ^ ] (3^) 



(3b) 



In Fig. 1 we show the time evolution of the CM mo- with t^fe = 2| sin(fc/2)|. We expand Eq. @ in Xn as: 



Xk = -ujlxk + - 'r— V [e''^'''^Xk,■■■Xk^Sk^+■■■+k„,.-q+k + {~l^e'''^'''^Xk^■■■Xk„Jk^+■■■+k„,,q+k\ (4) 
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Since in the ground state the only modes present in 
order A are Xq = x_q = A/2a;| the CM is coupled only 
to these modes up to second order in A: 

Q = iX-K (e'^'^^a;-, - e-^^^^a;,) (5a) 
Xq = -Lolxq + ^e^2-Q (5b) 

X.q = -Lulx.q + ^e-^2.Q (5,) 

In Fig. 2 we compare the behavior of Pit) — Q{t), ob- 
tained by solving the minimal set of Eqs. (H) with the ap- 
propriate initial conditions Q{t ~ Q) — 0, P{t = 0) = Pq, 
Xq{t = 0) = X/{2lj^), Xq{t = 0) = 0, with the one ob- 
tained from the full system of Eqs. (||). Eqs. (||) repro- 
duce very well the initial behavior of the CM velocity 
which displays oscillations of frequency A around the 
value r2/27r but do not predict the decay occurring at 
later times because, as we show next, this is due to cou- 
pling to other modes. To this aim, we analyze the relation 
between the initial CM velocity Pq and il/27r, respec- 
tively A±. 
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FIG. 2. Simulation of the full FK system according to Eq. 
(^) (solid lines) and numerical solution of Eqs. (^) (dashed 
lines) for TV = 144, A = 0.015, and several values of Pq around 
Po = 0.29. The differences between the two approaches 
are negligible. The average value of the CM velocity Q./2ti 
(horizontal dashed line) and the period of the oscillation for 
Po =0.29 are also shown. 
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Take as an ansatz for the CM motion: 



Q{t) 



2ti 



t + a+ sin(A+t) + sin(A_i) 



(6) 



Inserting the ansatz (||) in the coupled set of Eqs. (|^) 
keeping only terms linear in a±, we find that both 
are roots of: 



= A^TT^ (2Z(0) - Z(A) - Z(-A)) 
Z(A) being the impedance 
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In general Eq. has (besides the trivial solution A = 0) 
indeed two solutions, related to the sum and difference 
of the two basic frequencies in the system, ujq and f7 : 



A+ = 



(9) 



Close to resonance, VI ^ ujq, the amplitude a_ dominates 
(see below) and the CM oscillates with a single frequency 
A = A_ (see Fig. 2). Very close to resonance (more pre- 
cisely <Q. < Wq-|-(2A^7r^/wq) 3 ), the root A_ becomes 
imaginary, signaling an instability. In fact the system 
turns out to be bistable as it can be seen in Fig. 3 by the 
jump in r2(/o) ^ Pa passes through cjg/(27r). 




FIG. 3. Closed and open dots, frequencies A (left axis) 
and Q, (right a^cis) versus Po for simulations for A*" — 144, 
A = 0.015. Dashed and solid lines: solutions of Eqs. (^) and 
(^l|) respectively. 

Analytically, the relation between and Pq and the 
amplitudes a±, is determined by matching the ansatz 
(^ with the initial condition: 
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2^2 (cOq±n)^ 
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Pn 
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(n+UJq)^ {n-LJq)^ 



(11) 



The fact that Eq. ( |ll| ) has multiple solutions for 
when Pq ~ LOqj^-K is in accordance with the jump seen in 
Fig. 3. However, Eq. (pd|), which is derived by keeping 
only linear terms in Q!_, is not accurate enough to de- 
scribe in detail the instability in the above range around 
ujq where a_ diverges. 

An initial behavior similar to that for Pq — i^q/'i'K is 
observed in Fig. 1 for uPq ~ w„q/n27r. We examine the 
case n ~ 2. Eq. (^ shows that X2q is driven in next 
order in A by a;^: 



X2q = -UJ2qX2q + « 



«2A7re'2^'5a 



(12) 



When 2itQ ~ fit, Xq will be ~ Ae , so that X2q is forced 
with amplitude A^ and frequency 2$!, yielding resonance 
for 2ri = uj2q- Since X2q couples back to x^, we have a 
set of equations similar to Eqs. (||), but at order A^. 

We now come to the key issue, namely the onset of 
friction causing the decay of the CM velocity seen in Fig. 
1 at later times, which cannot be explained by the cou- 
pling of the CM to the main harmonics nq. Since Xq is 
by far the largest mode in the early stage, we consider 
second order terms involving Xq in Eq. (^: 
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q)] XU (13) 



Insertion of the solution obtained above for Xq (Eq. (ph|)) 
and Q (Eq. (|)) yields 



Xk = — [t^fe + A + B cos (At)] Xk 



(14) 



with A = 2(A7r)VZ(0) and B ^ a-. Clearly, Eq. (|l|) is 
a Mathieu parametric resonance for mode Xk- The rele- 
vance of parametric resonances has been recently stressed 
|]l2| . However here, due to the coupling of the CM to the 
modulation mode q, resonances are not with the wash- 
board frequency but with A ^ Q ~ ujq. Hence, we find 
instability windows around u'^+A — {nA/2)'^. Since A is 
small close to resonance, one expects to find instabilities 
for acoustic modes with k small. Indeed, as shown in Fig. 
4a, we find by solving Eq. (||) that the decay of the CM 
is accompanied by the exponential increase of the modes 
A: = 2, 3, 4 and, with a longer rise time, fc = 1. However, 
the instability windows resulting from Eq. ( [T^ ) , shown in 
Fig. 4b, cannot explain the numerical results of Fig. 4a, 
i.e. the Mathieu formalism cannot explain the observed 
instability. In Eq. (^, the only linear terms left out in 
Eq. ( p^ ) are couplings with Xk±q , which are much higher 
order in A. Nevertheless, these terms are crucial since 
they may cause new instabilities due to the fact that, for 
k small, they are also close to resonance. We have solved 
the coupled set of equations for mode x±k and Xk±q : 
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Xk = 
Xk±q = 



(15a) 
(15b) 



together with Eqs. (||) for continuous k. Indeed, we 
find a wider range of instabilities, giving a detailed ac- 
count of the numerical result as shown in Fig. 4b. This 
mechanism where a parametric resonance is enhanced by 
coupling to near resonant modes is quite general in sys- 
tems with a quasi continuous spectrum of excitations and 
is related to the one proposed in explaining instabil- 
ities in the FPU chain in a different physical context. 

The number of particles in the chain is an important 
parameter. When this number is very small, the chain is 
in fact commensurate and the phase of the CM is locked 
(the gap scales as due to Umklapp terms). Next, one 
enters a stage of apparent superlubric behavior due to the 
fact that the spectrum is still discrete on the scale of the 
size of the instability windows discussed above. For N — 
144 and A = ^Ac (Fig. |l|) we only begin to see the decay 
for values of Pq close to resonances. The experimentally 
observed superlubricity in fl^ could then be due either to 
the finiteness of the system or to the low sliding velocities. 

The above described multiple parametric excitation 
gives rise to an effective damping for the system via 
a cascade of couplings to more and more modes via 
the non-linear terms in Eq. (^). It remains an open 
question if this mechanism will eventually lead to a 
full or partial equilibrium distribution of energy over 
the normal modes JTot although our preliminary results 
support the former hypothesis even at weak couplings. 

In summary, we have described in detail the mecha- 
nism which gives rise to friction during the sliding of 
a harmonic system onto an incommensurate substrate. 
The onset of friction occurs in two steps: the resonant 
coupling of the CM to modes with wavevector related 
to the substrate modulation leads to long wavelength os- 
cillations which in turn drive a complex parametric reso- 
nance involving several resonant modes. This mechanism 
is robust in that it leads to wide instability windows and 
represents a quite general mechanism for the onset of 
energy transfer in systems with a quasi continuous spec- 
trum of excitations. 
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FIG. 4. (a) \xk{t)f of the first 4 modes from Eq. (|) with 
N = 144, A = 0.015 and Po = 0.29. Note that the first 
mode has a longer rise time and that the third mode is the 
most unstable, (b) Dispersion relation for a chain of 144 
atoms (fc-values in units of (^)). Unstable modes result- 
ing from the full simulation are represented by solid dots. 
The shaded fc-ranges give the instability windows resulting 
from the Mathieu-type Eq. (|l^) and cannot explain the simu- 
lation. Conversely the wiggled ranges of the phonon disper- 
sion (WU=weakly unstable, SU^strongly unstable) are the 
instability windows predicted by Eqs. (^) and (|l^). They ex- 
plain all instabilities as well as the long rise time of the first 
mode (WU) and the shortest one of the third mode which 
falls in the middle of the SU range. 
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